Antiferromagnetism and hole pair checkerboard in the vortex state of high Tc 

superconductors 
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We propose a microscopic state for the vortex phase of BSCO superconductors. Around the vortex 
core or above H C 2, the d wave hole pairs form a checkerboard localized in the antiferromagnetic 
background. We discuss this theory in connection with recent STM experiments. 
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Since the theoretical prediction of the antiferromag- 
netic (AF) order in the vortex state of the high T c super- 
conductors (SC)[1, 2], tremendous experimental[3, 4, 5, 
6, 7, 8, 9] and theoretical[10, 11, 12, 13, 14, 15, 16, 17, 
18, 19] progress has been made in recent years. Neutron 
scattering experiments[3, 4, 8, 20] in the Lai- x Sr x GuO± 
compound first reported enhanced magnetic scattering 
induced by the magnetic field in the vortex state. The 
field induced moment scales linearly (with possible loga- 
rithmic corrections) with the external field, or the num- 
ber of vortices in the system, in accordance with theo- 
retical predictions [2, 13]. The ^(SR experiment [5] in the 
underdoped YBa 2 Cu 3 Oj + s compound reveals static AF 
order in the vortex core. The NMR experiment. [6, 9] in 
the YBCO compound under high magnetic field shows 
enhanced AF ordering, predominantly from the region in- 
side the vortex core. These experiments have important 
implications on our theoretical understanding of high T c 
superconductivity. In the 5*0(5) theory, the transition 
from the AF Mott insulator phase to the SC phase is 
understood in terms of the rotation of the SO (5) super- 
spin vector as the chemical potential or doping level is 
varied[l]. Recent experiments in the AF vortex state 
show that the SO (5) rotation picture can be validated in 
real space, as a function of the distance from the center of 
the vortex core, where the superspin rotates continuously 
from AF to SC directions. 

While the salient features of the recent AF vortex ex- 
periments are in overall agreement with the theoretical 
predictions, there are additional experimental features 
worth exploring. A striking feature is revealed in the 
STM experiments by Hoffman et al.[7], where the local 
density of states (DOS) near the vortex core show a two 
dimensional checkerboard-like modulation with a 4a x 4a 
charge unit cell. Here a is the dimension of the CuO^ unit 
cell. This modulation decays exponentially away from the 
center of the vortex core, with a decay length of about 
10a. More recently, a similar pattern has also been seen 
in the absence of the applied magnetic field [21], possibly 
induced by the impurities at the surface. Understand- 
ing the four unit cell checkerboard modulation structure 
around the vortex core is the main focus of our present 



work. 

One possible explanation of the STM experiments is 
the formation of the stripe AF order [22, 23] around the 
vortex core or impurities. If the stripes are static, the 
question is how the ID stripes on two different layers 
of the bi-layer Bi-^Sr^C 'aC^Og+a material are oriented 
relative to each other. If we assume a strong coupling 
between the bi-layers, one would naturally expect the 
stripes on two different layers to be oriented along the 
same direction. This state breaks the rotational symme- 
try and is inconsistent with the 2D checkerboard pattern 
observed in the STM experiments. On the other hand, 
it is possible that the stripes on two different layers are 
oriented along the orthogonal directions. In this case, 
the STM tunneling should still predominantly probe the 
upper layer, and observe a strongly anisotropic quasi- 
1D pattern. Therefore, the static ID stripe picture can 
be experimentally falsified as the instrumental resolution 
improves. Current observations could be consistent with 
a picture of dynamically fluctuation stripes. The ques- 
tion is the time scale of the fluctuation. As the mag- 
netic field increases, these fluctuations are expected to 
slow down, and eventually, a static ID ordering pattern 
should emerge. This scenario could be tested in the un- 
derdoped materials or in the Nd2- x Ce x CuOi^s mate- 
rials, where H C 2 is within experimental reach. Besides 
the stripe state, one could also try to construct a Wigner 
crystal state of the doped holes in a AF background. 
Such an insulating state can only be realized for rational 
doping fractions. The optimally doped system is close to 
x = 1/8 doping, with one hole per eight sites. In this 
case, the holes can form a y/8a x \/8a superlattice uni- 
formly distributed over the AF spin background. Such 
a state has the full four-fold rotational symmetry as ob- 
served in the STM experiment. Unfortunately, the charge 
unit cell in this case is smaller than the experimentally 
observed 4a x 4a structure. Finally, the simplest expla- 
nation of the checkerboard charge structure around the 
impurities could be the Friedel oscillation associated with 
the fermi surface nesting [24] . The observed 4a x 4a charge 
pattern could be produced by the quasi-particles at the 
anti-node. However, these quasi-particles are gapped, 
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and it is not clear if they can explain the low energy 
feature at 7meV observed in the vortex core experiment. 

In this work, we propose a specific microscopic state 
around the vortex core and above H C 2 for the YBCO and 
BSCO class of superconductors. The insulating version 
of this state at doping x = 1/8 is depicted in Fig. 1 . In 
this state, the d wave hole pairs form a Wigner crystal, 
dubbed a d wave hole pair checkerboard. The center-of- 
mass of the d wave hole pair is located at the center of 
every four plaquette on the original lattice. This state 
has the 4a x 4a checkerboard symmetry as observed in 
the experiment, and the doping level for the insulating 
x = 1/8 state is reasonably close to the optimal doping 
level of the cuprates. Therefore, by forming the Wigner 
crystal state of the hole pairs rather than the holes them- 
selves, the doping level can be compatible with the ob- 
served size of the charge unit cell. Here the d wave pair 
is drawn on the plaquette for illustrative purposes. In 
the actual system, the pair size can be more extended. 
In this Letter, we discuss the case where the AF regions 
are connected, therefore, antiphase domain walls would 
be energetically unfavorable. However, the antiphase do- 
main walls would be possible if the hole pairs are more 
extended along the four-fold symmetry direction, and if 
there is an effective coupling among the triplet magnons 
on ncxt-nearest-neighbor plaquettes. We propose that 
the hole pair checkerboard state is energetically compet- 
itive in the microscopic t — J model with short ranged 
Coulomb interaction. This state takes advantage of the 
underlying AF exchange energy in the background, while 
the holes form singlet d wave bonds to minimize their ex- 
change energy as well. Finally, the kinetic energy of the 
hole pair competes with the Coulomb interaction. At 
x = 1/8, if the Coulomb interaction is stronger than the 
kinetic energy of the hole pair, the hole pair checkerboard 
state could win. We believe that the YBCO and BSCO 
materials are in the regime where the kinetic energy of 
the hole pairs favor the uniform SC state. However, when 
the applied magnetic field impedes the coherent hole pair 
motion either around the vortex core or in the state above 
H C 2, we argue that the hole pair checkerboard state can 
be realized, either in the static or dynamically fluctuat- 
ing form. Therefore, the central idea behind this work is 
that the magnetic field destroys the phase coherence of 
the hole pair by localizing them into a crystal, without 
breaking the pair. 

The hole pair checkerboard state arises naturally from 
the projected 5*0(5) model on a lattice[25, 26, 27]. In 
this approach, the original lattice is divided into non- 
overlapping plaquettes, depicted by solid squares in Fig. 
1. In the spirit of real-space renormalization, we retain 
only five low energy states on the plaquette, the half- 
filled singlet state, denoted by the triplet magnon 
state and d wave hole pair state t\\tt) with 2 holes 

on the plaquette. The d wave hole pair on a plaquette 
can be justified from the microscopic models [28], but a 



more extended hole pair can also be incorporated into 
our model. We define the center of the non-overlapping 
plaquettes to be the lattice sites of the projected SO (5) 
model, which has the spacing of 2a. Here we follow the 
notations and the conventions of reference [25] . The pro- 
jected 5*0(5) model is defined in terms of the following 
Hamiltonian. 

H = A s J2ttt a (x) + A c J2t{t h (x) 

X X 

{xx') 

- ^'£e i +<H) l (x)t h &) + h.c. 

{xx') 

+ V 1 J2 P(x)p(x') + V 2 J2 P( X M X ')> (!) 

(xx') ((xx')) 

where p(x) = t\th{x) is defined as the hole-pair density. 

As described in reference [25], the A s term describes 
the effective spin gap energy for the magnons, the A c 
term describe the effective chemical potential for the hole 
pairs, while the J s and J c terms describe the effective AF 
and SC Joscphson couplings respectively. Here e z ^ ij cor- 
responds to the phase introduced by the magnetic field 
in the vortex state, and we have also included the near- 
est neighbor and next nearest neighbor Coulomb inter- 
actions V\ and Vj2- This model can be approximately de- 
rived from the microscopic Hubbard and t — J model by 
a renormalization group transformation[26]. The model 
without the Coulomb terms have been studied exten- 
sively by numerical calculations, and the phase diagram 
has been deter mined [2 7]. 

If one neglects the spin degree of freedom, the 50(5) 
model reduces to the hard core boson model on a two 
dimensional square lattice. Even this simple model has a 
rich phase diagram as the density, the kinetic and the po- 
tential energies are varied[29, 30]. Besides the superfluid 
ground state, the system can also develop charge-density- 
wave order and supersolid order in the ground state. In 
the strong coupling limit at half-filling, the ground state 
develops checker-board charge order where every other 
site is occupied by a boson. At quarter-filling, a similar 
checker-board pattern develops where every four sites is 
occupied by a boson. This state is stabilized by the near- 
est neighbor and the next-nearest neighbor interactions 
V\ and V2. Including the AF spin order and translating 
the hole density and lattice size from the 50(5) lattice 
model to the original lattice, we find that the 1/4 filled 
checkerboard charge ordering of the hard-core bosons cor- 
responds exactly to the d wave pair checkerboard state 
illustrated in Fig. 1. However, we shall see that the AF 
exchange energies also play an important role in deter- 
mining the competition between the various states. 

We consider the following mean field state for the 
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FIG. 1: Illustration of the d wave hole pair checkerboard 
state at x — 1/8. In this state, the d wave hole pairs occupy 
every four non-overlapping plaquettes on the original lattice. 
The charge unit cell is 4a x 4a. The 5*0(5) model is defined 
on the center of the non-overlapping plaquettes. Such a state 
could be realized around the vortex core, whose center is de- 
picted by the cross. In the actual realization of this state, the 
hole pair can be much more extended, and the AF ordering 
could be much reduced from the classical value. 



Hamiltonian (1), 
I*) =Y[(cos6 + sin 9 (cos atl+ sin ae ill> tl))(x)\n} (2) 

X 

and minimize the energy E = as a function 

of the variational parameters 6(x),a(x) and 4>{x). At 
quarter filling of the hole pair bosons, (which corresponds 
to x = 1/8 in the actual system), the mean field theory 
predicts a quantum phase transition at V1+V2 = V c = 6 J 
for J c = 2J S = 2 J. For V1+V2 <V C , the system is in an 
uniform SC state. On the other hand, for V\ + V2 > V c , a 
hole pair checkerboard state depicted in Fig. 1 is favored. 
We shall choose our parameters such that the system is 
in an uniform SC state, but it is close to the hole pair 
checkerboard state. 

In order to reduce the number of variables, we have 
assumed four-fold rotational symmetry as well as sym- 
metry by reflection along the diagonals. We can further 
simplify the problem by fixing the SC phase <f> assuming 
that the magnetic field is uniform in a region of radius 
equal to the London penetration depth (here, we present 
data with 30a but we have checked that this effect is very 
small). 

The boundary sites are taken to be in a bulk SC phase 
with a fixed hole-pair density, which is taken to be the 
same as the averaged one. We fixed the center of the 
vortex core to be fully anti- ferromagnetic, i.e., a = 
and 9 = ir/4. 

In Fig. 2, we present the result of charge distribution 
which shows a clear checkerboard pattern. The contrast 
decreases very quickly away from the vortex core over 



a distance ~ 10a which is roughly consistent with the 
experimental value [7] . We would like to emphasize that 
this checkerboard pattern persists for smaller interactions 
V\ and V2, but with a smaller contrast. For sufficiently 
small contrast, the modulation of the charge density may 
be hard to detect in the integrated LDOS[7]. 
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FIG. 2: Charge distribution around the vortex core. Pa- 
rameters of the Hamiltonian are: Vi/J = 3, V2/J = 2.2, 
A s /J = 1.2 and A c /J = -0.12. The center of the figure is 
the center of the vortex. The grid size is 2a. 

We also present the profiles of the AF order parame- 
ter (n 2 ), the SC order parameter (ni) and the 7r order 
parameter (iL 2 5)/2[l] in Fig. 3. Notice that the w order 
parameter is only non-zero in the region where the AF 
and SC order parameters overlap. In our projected SO (5) 
model, the 5*0(5) orthogonality condition]!, 14] 



(L a b)(n c ) + (L bc )(n a ) + (L ca )(n b ) = 



(3) 



is satisfied at every site. This relation predicts a quanti- 
tative relation among these three order parameters and 
the hole density, p(n 2 ) + (£25) = 0, which is satisfied 
by our numerical results of the order parameters. The 
relative size of the AF and SC regions is determined by 
the ratio of J s /J c , and can be varied to account for the 
experimental observations at low temperature. 

We shall now discuss the experimental consequences 
of our theoretical proposal. Our proposal can be tested 
most directly in the vortex state, or above H C 2, by mea- 
suring the static charge structure factor corresponding to 
a 4a x 4a unit cell with XY rotational symmetry. The 
magnetic structure factor could either be a commensu- 
rate (tt, 7t) peak with an amplitude modulation of four 
unit cells or a shifted peak corresponding to an anti-phase 
modulation. Since the energy difference between these 
two types of magnetic order is very small, we conjecture 
that both types could occur in difference compounds of 
the high T c cuprates. Our proposal differs in detail from 
the stripe picture from the fact that the ordering pat- 
tern has square symmetry and from the fact that the 
ordering charge is the hole pair rather than single holes. 
Our proposed state is also consistent with the ^.sR and 
the NMR experiments in the vortex state of the YBCO 
system[5, 9]. The fxsR experiment reports commensu- 
rate (ir, 7r) AF order inside the vortex core[5]. Because 
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FIG. 3: From top to bottom: Distribution of the SC, AF and 
7r order parameters around the vortex core. The relative size 
of the AF and SC regions can be varied by varying Js/Jc- 
Notice that the n order parameter is only non-zero in the 
region where the AF and SC order parameters overlap. 



of the modulation of the AF order parameter amplitude, 
oxygen sites, for examples those located in between the 
singlet hole pair and the AF plaquettes, can still detect a 
spin signal [5, 9]. In this work we primarily focused on the 
charge ordering in the vortex state. Similar arguments 
may be applied to the impurity states as well [21]. How- 
ever, we caution here that if the impurity is pair breaking, 
the idea of a order state of hole pairs may not be directly 
applicable. We would also like to call attention to ref. 
[31], where experimental evidence for Wigner lattice was 
reported. 

We also encourage the detailed study of this hole 
pair checkerboard state in the numerical calculations 
of the t — J model, including effects of the Coulomb 
interaction[32]. While we do not propose this state to 
be realized in the ground state for physical range of pa- 
rameters in the absence of the magnetic field, it is inter- 
esting to test if this state is realized in the vortex state 



or above H C 2- It is also instructive to determine param- 
eters for which this state is realized in the ground state, 
and see how far these parameters are from the physically 
reasonable values. This way, we can estimate the pertur- 
bations required for this state to reveal itself. 
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